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Abstract: In nature, the mechanical properties of geological bodies are very complex, and its 
various mechanical parameters are vague, incomplete, imprecise, and indeterminate. In these cases, 
we cannot always compute or provide exact/crisp values for the joint roughness coefficient (JRC), 
which is a quite crucial parameter for determining the shear strength in rock mechanics, but we need to 
approximate them. Hence, we need to investigate the anisotropy and scale effect of indeterminate JRC 
values by neutrosophic number (NN) functions, because the NN is composed of its determinate part 
and the indeterminate part and is very suitable for the expression of JRC data with determinate and/or 
indeterminate information. In this study, the lower limit of JRC data is chosen as the determinate 
information, and the difference between the lower and upper limits is chosen as the indeterminate 
information. In this case, the NN functions of the anisotropic ellipse and logarithmic equation of JRC 
are developed to reflect the anisotropy and scale effect of JRC values. Additionally, the NN parameter 
p is defined to quantify the anisotropy of JRC values. Then, a two-variable NN function is introduced 
based on the factors of both the sample size and measurement orientation. Further, the changing 
rates in various sample sizes and/or measurement orientations are investigated by their derivative 
and partial derivative NN functions. However, an actual case study shows that the proposed NN 
functions are effective and reasonable in the expression and analysis of the indeterminate values 
of JRC. Obviously, NN functions provide a new, effective way for passing from the classical crisp 
expression and analyses to the neutrosophic ones. 


Keywords: joint roughness coefficient (JRC); anisotropy; scale effect; neutrosophic number; 
neutrosophic number function 





1. Introduction 


According to the shear strength formula given by Barton [1] in 1973, the joint roughness 
coefficient (JRC) obtained effectively is quite a crucial parameter for determining the shear strength 
in rock mechanics. Various methods for determining the JRC values have been proposed by many 
scholars [2-10] since Barton [11] firstly introduced the concept of JRC. Since then there have also been 
many studies about the anisotropy and scale effect, which are two important characteristics of JRC 
in the sample sizes and measuring orientations. Barton and Bandis [12] proposed that the scale effect 
of roughness coefficient is one of the main reasons for the scale effect of mechanical behavior of 
rock mass through model test studies. Du et al. [13] found that the undulation amplitudes implied 
the negative scale effect and presented a good logarithmic function relationship with alternation 
sizes. Recently, Chen et al. [14] used the digital image processing (DIP) technique to retrieve the joint 
surface topography and a variogram function for characterizing the anisotropy of the joint surface 
roughness and to estimate the joint aperture. However, existing research methods can only deal 
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with the determinate/crisp function expressions of JRC values and analyze their anisotropy and size 
effect, but cannot express and handle the incomplete, uncertain, and imprecise data of JRC. Due to 
the incompleteness of our observations, measurements, and estimations, or the existing disturbances 
and uncertainties in the real world, we cannot always compute or provide exact/crisp values for 
the inhomogeneity of the rock joints, but we need to approximate them in indeterminate environments. 
In this case, Ye et al. [15] established neutrosophic functions (interval functions) of JRC and shear 
strength. However, they only gave the interval/thick functions of JRC and shear strength, but cannot 
express such an indeterminate function containing coefficients and/or parameters of neutrosophic 
numbers (NNs). 

In fact, there always exists determinate and indeterminate information in the real world. Hence, 
Smarandache [16-18] firstly introduced the concept of NN, which is really suitable for expressing 
determinate and/or indeterminate information, denoted as z = a + bI. Obviously, the NN z consists 
of its determinate part a and its indeterminate part bI, where a and b are real numbers and I denotes 
the indeterminacy. Clearly, NN is more suitable for the expression of determinate and/or indeterminate 
information in the real world. In recent years, Ye [19] firstly proposed a multiple attribute group 
decision-making method under a NN environment and developed the possibility degree ranking 
method for ranking alternatives. Then, Ye [20] presented a bidirectional projection model of NNs for 
multiple attribute group decision-making problems with NN information. Further, Kong et al. [21] and 
Ye [22] proposed the cosine similarity measure of NNs for the misfire fault diagnosis of gasoline engines 
and the exponential similarity measure of NNs for the fault diagnosis of steam turbines under a NN 
environment. Additionally, Chen and Ye [23] introduced a projection model of NNs for the multiple 
attribute decision-making problem of clay-brick selection under a NN environment. Further, Liu and 
Liu [24] presented a NN generalized weighted power averaging operator and applied it to multiple 
attribute group decision-making under a NN environment. Ye [25] presented a NN linear programming 
method for NN linear programming problems under an indeterminate environment. All of these 
studies have shown the superiority of NNs to express indeterminate information. Unfortunately, 
NN functions have not been studied and applied in rock mechanics thus far. Hence, this original 
study will establish NN operations and NN functions of JRC with determinate and/or indeterminate 
information to express and analyze the anisotropy and scale effect of indeterminate JRC values and 
their changing rates in measuring sizes and orientations of test samples by the actual case of the natural 
joint surface sample obtained from Changshan County of Zhejiang province in China. In the real 
world, we cannot always compute or provide determinate/crisp JRC values for the JRC characteristics, 
but we can approximate them to realize the indeterminate expression and analysis of JRC values. 
In this case, the NN functions provide an effective mathematical tool for passing from the classical 
crisp expression and analysis to the neutrosophic ones. 

The rest of this paper is organized as follows: Section 2 presents some basic concepts and 
operations of NNs; Section 3 introduces the JRC data acquisition of natural rock joint surfaces 
with an actual case; In Section 4, some elliptic equations with NNs (elliptic NN functions of JRC) 
are introduced for the descriptions of the joint profiles in various orientations and their changing 
rates are discussed by their derivative functions, and then the NN parameter y is defined to quantify 
the degree of the anisotropy; In Section 5, logarithmic equations containing NNs (logarithmic NN 
functions of JRC) are presented to analyze the JRC values and the changing rates in differently-sized 
joint samples; Section 6 introduces a two-variable NN function of JRC to describe JRC values in various 
orientations and sizes and discusses its changing rates by its partial derivative functions; Section 7 gives 
discussion; Finally, conclusions and future work are given in Section 8. 


2. Some Basic Concepts and Operations of NNs 


The concept of NN was firstly proposed by Smarandache [16-18], which consists of a determinate 
part and an indeterminate part. Then, it is denoted as z = a + bI, where a and b are real numbers and I 
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is the indeterminacy, implying some indeterminate range I e [I‘, IY]. It is obvious that it can describe 
determinate and/or indeterminate information. 

For example, assume that a NN is z = 8 + 5I, where I e [0, 0.5]. Thus, its determinate part is 8, 
its indeterminate part is 5I, and then there is z e [8, 10.5] for I e [0, 0.5]. Furthermore, one can also 
specify some indeterminate range I e [I/, IY] to satisfy the actual applied requirements. 

Let zı =a, + bil and z2 = az + bzI be two NNs for I e [I, IŻ]; their operational relations are given 
as follows [25]: 





























zı +22 = ay +a + (b1 + bo) = [ay +02 + bI" + bol”, ay +a + bI" + byl] (1) 
z1 — 22 = ay — m + (by — by) I = [ay — az + b1 I% — Dg", ay — ay + b1 IY — b21] (2) 
Z1 XZ? = ma + (aby + azb1)I + b1b21? 
min (a, +b1I}) (az + bol"), (ay + by") (az + bel”), 
_ (a1 + by 1) (az + bol”), (ay + b11Y) (ag +b:1€) j’ (3) 
ae (ay + by 1") (ap + bol"), (a1 + bI!) (az + byI¥), 
(a, + by I€) (a + bo”), (a1 + by I") (ag + bg I") 
zy ay thy [ay tb Thay +b, 14] 
22 a2+b21 [ay +b2I},a2+b21U] 
. ay +b, I} ay +b, IL ay +b, 1U ay +b, IU 
min az+b21U” ay +b Ih’ ay+by IU’ ay +bo IE d (4) 
= max( a, +b, I} ay +b, IE ay +b, IU att) 
a+b IU’ a7 +b TE” a+b IU’ a7 +by TE 











3. JRC Data Obtained from an Actual Case 


In real situations, we choose a natural joint surface sample to collect data from Changshan 
County of Zhejiang province in China as an actual case. The foliated slate rock is grayish, completed, 
and hard and has a dense structure. Its surface undulates from 1 cm to 2 cm and is lightly weathered. 
Then, a simple mechanical hand profilograph invented by Du [26] was employed to record the surface 
profiles by every interval of 15° in [0°, 360°] along the specified lines. Then, we intercepted samples 
from 10 cm to 100 cm at an interval of 10 cm in every orientation. Finally, we use the models proposed 
by Zhang et al. [27] to calculate the JRC values of each samples in different sizes and directions. 
By the statistical analysis of all JRC values of every sample, we can obtain the average values and 
standard deviations of all sizes in each direction, which are used in the following case analyses. 


4. NN Functions for Describing the Anisotropic Characteristics of JRC 


4.1. NN Functions of JRC in All Orientations 


Here, we take a sample in the length of 20 cm as an example. All measured values and the mean 
values of JRC in every orientation on one polar plot are shown in Figure 1. It can be found that 
the JRC values are randomly distributed from 0 to 20 in each orientation. Then, the JRC values 
in the orientations of 90° and 270° are smaller, while the JRC values in some orientations such as 0° 
and 180° are relatively greater. Performing a calculated analysis, the JRC average value and variance of 
each angle orientation are obtained by the statistical analysis. For example, the average value JRCimean 
and standard deviation ø of JRC values for 20 cm joint samples in the 0° orientation are 9.653 and 1.716, 
respectively. The smaller standard deviation indicates that the data points tend to the mean value. 

According to the statistical analysis, about 70% of the points focus on one area in Figure 1, then most 
of the blue stars are crowded together. Hence, we may choose a robust range (confidence range) from 
JRCmean — 7(lower limit) to JRCmean + o(upper limit) for describing the anisotropic characteristics. 
In Figure 1, it is shown that these upper and lower limits can be approximately structured as one ellipse 
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area. It is obvious that the JRC mean values and most of the sample data fall into the area between two 
red ellipses. 

Based on the least square method, the upper and lower limits of JRC values are fitted by two 
ellipse envelopes for joint samples of each size. Then, according to the concept of NNs and the general 
parametric equations of the ellipse, the NN parametric equations of JRC can be written as: 


„I € [0,0.5] (5) 





x = (a, + bI) cosé 
y = (a2 + bol) sin 0 


where a, and a are two radiuses of the inside ellipse envelope, and a, + bI and az + bol are two 
radiuses of the outside ellipse envelope. Through the statistical analysis of all sizes, the indeterminacy 
I belongs to the interval [0, 0.5]. Then, b; (i = 1, 2) can be determined by the twice of the difference 
between the radii of the inside and outside ellipse envelopes. 

Each JRC value in Figure 1 is the distance from any point within the two ellipses to the coordinate 
origin. To express the JRC values in any orientation, by using Equations (3) and (5) we can derive 
the following NN function: 





JRC(0) = fx? +y2 = Vila + b11) cos 0}? + [(a2 + b21) sin 0}? 

= „/ (a1 + b11) + (a1 + bil) cos? 0 + (a2 + b21) - (a2 + bal) sin? 0 

= „/ a2, (a1 + b1/2)?] cos? 6 + [a3, (az + b21/2)?] sin? 0 

= [ya cos? 0 + a2 sin? 6, via + b,/2)* cos? 0 + (az + 2/2) sin? 6]for I € [0,0.5]. 








(6) 








In order to reflect the changing rate of JRC(@) in any orientation; here, we carried on the derivative 
operation of Equation (6) and obtain the following derivative Equation (7) so as to analyze the variation 
situations of anisotropy: 


JRC' (0) = | (a5 = a?) cos sin @ ((az + byI)* — (ay £ byI)*) cos bain’ 


a2 cos? 0 + a3 sin? 0 Va + b11)? cos? 0 + (az + b21)? sin? 0 








(7) 


For example, when the length L = 20 cm in Figure 1, the NN parametric equations with I e [0, 0.5] 
can be obtained by the following steps: 


(1)a, = 8.3,a2 = 5.7, 
(2)b, = (11.5 — 8.3) x 2 = 6.4, b2 = (9 — 5.7) x 2 = 6.6, 
sf x = (8.3 + 6.41) cos 0 


for I 5]. 
y= 766d T Eaa] 





Thus, we have the following NN function with I e [0, 0.5]: 





JRC(0) = 2 +y? = y [(8.3 +6.41) cos 6]? + [(5.7 + 6.61) sin 0]? 

= (8.3 + 6.41) - (8.3 + 6.41) cos2 9 + (5.7 +6.61) (5.7 +6.61) sin? 6 

= 1/832, (8.3 + 6.4/2)?] cos? 6 + [5.72, (5.7 + 6.6/2)?°] sin? 0 

= [v8.32 cos? 0 + 5.72 sin? 0, (83 + 6.4/2)? cos? 6 + (5.7 + 6.6/2)* sin? 6] 
= [68.89 cos? 6 + 32.49 sin? 6, \/132.25 cos? 6 + 81 sin? 6]. 
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Thus, its NN derivative function with I e [0, 0.5] is given as follows: 


—36.40 cos 0 sin 0 —51.25 cos 0 sin 0 


JRC' (0) = ; 
2 + ad 
\/ (68.89 cos? 0 +32.49sin29 v 132.25 cos? 8 + 81 sin? 6 








The changing rate of JRC’(@) is shown in Figure 2. It is obvious that the changing curves of 
the inside and outside ellipse envelopes present periodicity and the changing rates fall into the range 
between two negative sine functions. 





JRC of L = 20 cm,orientation(°) 
test samples 
—$— JRC _ +0 90 45 
mean pane, 
—_@—JRC-o 120. ~ 60 
mean a $ h S . 
Outside ! : 
ellipse envelope : 
inside i 


ellipse envelope 
JRC values 
mean 





Figure 1. Polar plots of directional JRC values, upper and lower limits of JRC values and 
the corresponding ellipse fitting envelopes of 20 cm joint samples. 





3 L = 20 cm,the rate of change 


— Inside ellipse 
oL — Outside ellipse a 





= 











-3 I li li i li li li 
0 50 100 150 200 250 300 350 400 
A(°) 


Figure 2. The changing curves of the lower and upper limits for JRC’(6). 


By the similar method, we can obtain the NN functions of JRC(@) and JRC’(@) with I e [0, 0.5] and 
various scales. Then, all of the results are shown in Tables 1 and 2. 
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Table 1. NN functions of JRC(@) for the JRC expressions of different sized joint samples and I e [0, 0.5]. 















































Scale (cm) ay bı a2 bo JRC(0) 
10 86 78 58 84  [/73.96cos?@ +33.64sin? 0, 156.25 cos? 0 + 100 sin? 6] 
20 83 64 57 66 — [v68.89 cos? 0 + 32.49 sin? 0, vV132.25 cos? 0 + 81 sin? 6] 
30 8 56 55 6 — [/64cos?@ +30.25sin2 0, v116.64cos2? 0 + 72.25 sin? 6] 
40 7.7 56 53 54 — [/59.29cos? 6 + 28.09 sin? 0, vV110.25 cos? 0 + 64sin? 6] 
50 7.4 5.8 5.2 4.6 [V 54.74 cos? 6 + 27.04 sin? 0, V/106.09 cos? 0 + 56.25 sin? 6] 
60 73 58 5.3 5.0 [53.29 cos? 6 + 28.09 sin? 6, v/104.04 cos 0 + 60.84 sin? 6] 
70 72 58 5.3 46  [\/51.84cos? 6+ 28.09 sin? 6, V/102.01 cos? 0 + 57.76 sin? 6] 
80 70 56 52 44 — [/49cos? @+27.04sin2 0, V/69.04 cos? 0 + 54.76 sin? 0] 
90 7.0 5.0 5.1 3.8 [V49 cos? 6 + 26.01 sin? 0, \/90.25 cos? 0 + 49 sin? 6] 
100 70 44 51 34 — [y/49cos?@ + 26.01 sin? 0, V/84.64 cos? 0 + 46.24 sin? 6] 








Table 2. NN functions of JRC’(@) of different sized joint samples and I e [0, 0.5]. 













































































cale (cm 
Scale (cm) RC’(0) 
10 —40.32 cos @ sin 0 —56.25 cos ĝ sin 
\/ (73.96 cos? 6+33.64 sin? 9” v/156.25 cos? 0+100 sin? 8 
—36.40 cos @ sin @ —51.25 cos @sin@ 
j 
20 \/ (68.89 cos? 0+32.49 sin? 6 ” 132.25 cos? 0+81 sin? 0 
—33.75 cos @ sin @ —44.39 cos @ sin @ 
30 v (64 cos? 0+30.25sin? 9 ” v/116.64cos? 0+72.25 sin? 6 
—31.20 cos @ sin 0 —46.25 cos ð sin 0 
40 / (59.29 cos? 6+28.09 sin? 9 ” v/110.25 cos? 0+64sin? 8 
—27.70 cos @ sin @ —49.84 cos ð sin 
50 y (54.74 cos? 0+27.04 sin? 9 ” 106.09 cos? 0-+56.25 sin? 8 
—25.2 cos @ sin 0 —43.2 cos @ sin 0 
60 y (53.29 cos? 0+28.09 sin? 0 ” v/104.04 cos? 0+60.84 sin? 6 
—23.75 cos @ sin 0 —44.25 cos ĝ sin 
70 \/ (51.84 cos? 6+28.09 sin? @ ’ ¥//102.01 cos? 0+57.76 sin? 6 
—21.96 cos @ sin 0 —14.28 cos @ sin @ 
80 y (49 cos? 0+27.04sin? 9 ” v/'69.04 cos? 0+54.76 sin? 0 
—22.99 cos 0 sin 0 —41.25 cos ð sin 
90 (9 cos? 0+26.01 sin? 0  v/90.25 cos? 0+49 sin? 0 
—22.99 cos @ sin 0 —38.4 cos ĝ sin 0 
100 y (49 cos? 0+26.01 sin? 0 ’ /84.64 cos? 0-+46.24 sin? 0 


4.2. Anisotropic Parameter 


In 2009, Tatone [28] defined a parameter to quantify the anisotropy that is the ratio of 
the maximum directional roughness to the minimum directional roughness on the polar plots. Based on 
the anisotropic definition, we present the NN anisotropic parameter to characterize the degree of 
anisotropy for each size. Here, the maximum directional roughness is the semi-major axis of ellipse; 


the minimum directional roughness is the semi-minor axis of the ellipse. Hence, the NN anisotropic 


parameter is defined as the following form: 


AT bıl 
fa pn (0, 0.5] 
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Then, according to the Equation (4), the anisotropic parameter y is calculated by the following formula: 











. (ae ay+b,x0 a,+b,x0.5 ayt+b,x0.5 
ay+b 7 x0.57 ag+b2 x0’ ag+b2x0.5’ ag+b2x0 /” 
yp = ani = max( ay by x0 ai+bıx0 ay+b,x0.5 aes) 
awe az +b X0.57 ay +by XO" a+b X0.5" a+b x0 (8) 





= 2 2aı +b 
[y a; 1,1 = [0,0.5]. 


For example, when L = 20 cm, we have the following result: 





3+6. ; 3 +64 
83 s] 2x83 2x83+64] 10.92 2.02], I € [0,05] 


Y 5.7 + 6.61 2x5.7+66' 2x57 





Based on the NN anisotropic parameter y, the interval ranges of the anisotropic parameters of all 
sizes are shown in Figure 3. Both the upper and lower bounds are fitted by the logarithmic curves, 
where the upper curve decreases and the lower curve increases as L increases. 














2.2 
Mi 7 -Up 
2F z + + + Wdown - 
1.8- ż i A | 
W p = 2.478-0.136In(L), R = 0.9471 
1.6} 1 
S 
1.44 + 
1.2 w= 0.7427+0.05457In(L), R= 0.8392 , 7 
down 
1F % + z = ia = 7 
0.8 ii li li li li li li i li li 
0 10 20 30 40 #50 #60 70 #980 90 100 110 
L Cem) 


Figure 3. The upper and lower bounds of the parameter y for different sized joint samples. 


From Figure 3, we can obtain the NN anisotropic parameter (L) corresponding to L for I e [0, 0.5] 


as follows: 


(L) = 0.7427 + 2.4721 — (0.05457 + 0.162861) In L 0) 
= [0.7427 — 0.05457 In L, 2.478 —0.136In L], I = [0, 0.5] 


5. NN Functions for Reflecting the Scale Effect of Joint Surface Roughness 


As we know, the roughness coefficient of a rock mass structure decreases with the increase of size, 
which is called the scale effect. In this section, we discuss this issue. 

When the orientation is 15°, the JRC values on the increase of the joint sample sizes are shown 
in Figure 4, which apparently exhibit a negative scale effect. Then, the fitting curves present 
the logarithmic functions with monotonously-decreasing trends, and then both the upper and lower 
bounds are matched well. 

Thus, the NN function JRC(L) and its derivative function JRC’(L) for I e [0, 0.5] are given as: 


JRC(L) = (a, +b1I) + (a2 + b21) NL = (a, + a2 nL) + (bı +b21n L)I 


10 
= [ay + ag InL,a, + aln L+ (by + by In L)/2],1 E [0,0.5], g) 


bol 
JRC’(L) = Or, 


I € [0,0.5] (11) 
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where a, and a are the fitting parameters of the lower bounds and bı and b are the twice of 
the difference of the lower and upper bounds. In the NN function JRC(L), 41 + azInL is the determinate 
part and (b1 + b2lnL)I is the indeterminate part. 
For example, in Figure 4, a; and az are 9.375 and —0.3999, respectively. Then, bı and b} can 
be obtained as follows: 
by = 2 x (16.74 — 9.375) = 14.73 
by = 2 x (—1.653 + 0.3999) = —2.5062 


According to Equations (10) and (11), we have the following results: 


JRC(L) = 9.375 — 0.3999 In L + (14.73 — 2.5062 In L)I 
= [9.375 — 0.3999 In L, 16.74 — 1.653InL], @ = 15° 


_ 0.3999 — 2.50621 | 1.653 0.3999 
L L” L 








JRC' (L) 


In Figure 5, the shaded part is the changing rate range of JRC’ (L). It is monotonously increasing 
on the increase of L. Clearly, the scale effect phenomenon is more obvious in smaller sizes. 

By the similar method, fitting results of other orientations can also be obtained, the fitting parameters, 
NN functions for JRC(L), and derivative functions for JRC’(L) are tabulated in Tables 3 and 4. 


d= 15° 
14 T T T T 








A JRC eant? 
JRC p = 16.74-1.653*InL,R® = 0.9669 
p JRC 


mean 


JRC 


mean 
fitting curve (up) 
fitting curve (down) 











o 

















8- | 
= 9/375-0.3999"InL,R* |= 0.891 
















































































0 20 40 60 80 100 120 
L(cm) 


Figure 4. JRC values of different sized joint samples in the orientation of 0°. 
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JRC'(L) 
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“0 10 20 30 40 50 60 70 80 90 100 110 
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Figure 5. The changing rates of upper and lower limits of JRC’(L) in various sizes. 
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Table 3. NN functions of JRC(L) for I e [0, 0.5] corresponding to the JRC values of the joint samples 


in different orientations. 























Orientation (°) a4 by a bz JRC(L) 

0 9.487 11.246 —0.5415 —1.253 [9.487 — 0.5415 In L, 15.11 — 1.168 In L] 
15 9.375 14.73 —0.3999 —1.653 [9.375 — 0.3999 In L, 16.74 — 1.653 In L] 
30 9.22 15.34 —0.4101 —2.5498 9.22 — 0.4101 In L, 16.89 — 1.685 In L 
45 8.83 12.78 —0.4538 —1.8364 8.83 — 0.4538 In L, 15.22 — 1.372In L 
60 7.432 14.736 —0.3859 —2.1142 7.432 — 0.3859 In L, 14.8 — 1.443 In L 
75 6.454 13.952 —0.2913 —2.3814 6.454 — 0.2913 In L, 13.43 — 1.482 In L 
90 4.447 15.086 0.1112 —2.5844 4.447 + 0.1112 1n L, 11.99 — 1.181 In L 

105 6.18 13.62 —0.2751  —2.1678 6.18 — 0.2751 In L, 12.99 — 1.359 In L 
120 7.748 15.544 —0.359 —2.712 7.748 — 0.359 In L, 15.52 — 1.715 In L 
135 9.009 10.482 —0.6697 —1.4586 9.009 — 0.6697 In L, 14.25 — 1.399 In L 
150 9.454 15.032 —0.5928 —2.4184 9.454 — 0.5928 In L, 16.97 — 1.802 In L 
165 10.49 12.58 —0.823 —1.746 10.49 — 0.823 In L, 16.78 — 1.696 In L 
180 8.705 12.55 —0.4852 —1.6456 8.705 — 0.4852 In L, 14.98 — 1.308 In L 
195 8.905 14.79 —0.5113  —2.3994 8.905 — 0.5113 In L, 16.3 — 1.711 In L 
210 9.395 11.33 —0.7162 —1.5076 9.395 — 0.7162 In L, 15.06 — 1.47 In L 
225 8.614 9.932 —0.6189 —1.1822 8.614 — 0.6189 In L, 13.58 — 1.21 In L 
240 7.112 10.756 —0.4301 —1.6358 7.112 — 0.4301 In L, 12.49 — 1.248 In L 
255 5.94 11.88 —0.2837  —1.8705 5.94 — 0.2837 In L, 11.88 — 1.219 In L 
270 4.659 12.042 —0.0052  —1.7540 4.659 — 0.0052 In L, 10.68 — 0.882 In L 
285 6.779 6.502 —0.5106  —0.5216 6.779 — 0.5106 In L, 10.03 — 0.771 In L 
300 7.978 10.324 —0.5043 —1.5754 7.978 — 0.5043 In L, 13.14 — 1.292 In L 
315 9.2 14.78 —0.6081 —2.2038 [9.2 — 0.6081 In L, 16.59 — 1.71 In L] 

330 8.792 13.676 —0.5216 —2.1808 8.792 — 0.5216 In L, 15.63 — 1.612 In L 
345 9.491 12.583 —0.6661 —2.0618 9.491 — 0.6661 In L, 15.76 — 1.697 In L 








Table 4. NN functions of JRC’(L) for I e [0, 0.5] in different orientations. 


Orientation (°) 


JRC'(L) 





0 
15 
30 
45 
60 
75 
90 

105 
120 
135 
150 
165 
180 
195 
210 
225 
240 
255 
270 
285 
300 
315 
330 
345 


—1.168/ L, —0.5415/L 
—1.653/ L, —0.3999 / L 
—0.4101/L, —1.685/ L 
—1.372/ L, —0.4538/L 
—1.443/ L, —0.3659/L 
—1.482/ L, —0.2913/ L 
[—1.181/L,0.1112/L] 
—1.359/L, —0.2751/L 
—1.715/ L, —0.359/L 
—1.399/ L, —0.6697 / L 
—1.802/ L, —0.5982/L 
—1.696/ L, —0.823/L 
—1.308/ L, —0.4852/L 
—1.711/L, —0.5113/ L 
—1.47/ L, —0.7162/ L 
—1.21/L, —0.6189/ L 
—1.248/ L, —0.430/L 
[—1.219/L, —0.2337/L] 














[—0.8822/L, —0.005181/ L] 
[—0.7714/L, —0.5106/L] 


[—1.292/L, —0.5043/L] 
[-1.71/L, —0.6081/L] 
[—1.612/L, —0.5216/L] 
[—1.697/L, —0.666/ L] 
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6. NN Function of JRC with the Two Variables 0 and L 


As analyzed above, we have obtained the anisotropy and the scale effect of JRC, respectively, 
which are the important characteristics of the rock joint samples. In this section, we shall describe 
a NN function with the two variables @ and L, under NN environments. 

To find the relationship between the JRC values and the two variables 0 and L, through numerical 
analyses, the trigonometric function polynomial of 0 and L can fit into two smooth surfaces of the upper 
and lower bounds in the three-dimensional coordinate system, as shown in Figure 6. It is obvious that 
all of the points of mean values fall between the upper and lower surfaces. That is to say, the uncertain 
range almost contains all JRC values. Thus, the fitting equations based on the upper and lower surfaces 
of JRC values are presented as follows: 


JRC(0, L) gown = —0.4536InL + 1.382sin0 + 3.866sin2@ + 1.039cos*@ — 32.03sin*6 — 
22.65sin°6 + 85.86sin°@ + 53.45sin’6 — 109.4sin8@ — 36.4sin?@ + 60.38sin!°@ — 
7.551cos!!@ — 11.04sin!26 + 6.802cos!9@ + 4.334sin!50 + 8.869, 


JRC(0, L)up = —1.42InL + 2.08sin@ — 9.352sin20 + 2.061cos?9 + 71.78sin*6 — 46.69sin°@ 
— 255.3sin°@ + 117.2sin’6 + 444.2sin®@ — 81.64sin?@ — 383.1sin!°@ — 18.68cos!!6 + 
128.9sin!20 + 16.96cos!36 + 9.156sin!°6 + 15.74. 


Then, the NN function JRC(@, L) for I e [0, 0.5] is given as follows: 


JRC(6, L) = (—0.4536 — 1.9328])InL + (1.382 + 1.396])sin@ + (3.866 — 26.4361)sin26 + 
(1.039 + 2.044I)cos*@ + (—32.03 + 207.62I)sin*6 + (—22.65 — 48.08I)sin°@ + 
(85.86 — 682.32])sin°@ + (53.45 + 127.5])sin’6 + (—109.4 + 1107.2/)sin86 + 

(—36.4 — 90.48I)sin?9 + (60.38 — 886.96])sin!6 + (—7.551 — 22.258]I)cos!!6 + 
(—11.04 + 279.88])sin!*6 + (6.802 + 20.316I)cos!70 + (4.334 + 9.644])sin@ + 
(8.869 + 13.7421), for I = [0, 0.5]. 


(12) 


i IRC down 
a] JRC p 


mean 


100 









40 L(cm) 


Figure 6. Curved surface area of JRC values corresponding to both L and 0. 


By the partial derivative of Equation (12), we can obtain the changing rate functions of the scale 


effect and the anisotropy: 
AJRCup (OL) _ 
aL = 





1 
1.42Ł, 


13) 
IJRCaown (O,L 1 ( 
domn(@L) _ _ 9.45361, 





Information 2017, 8, 69 11 of 13 


PIRC (eL) L 9.08 cos 6 — 18.704 cos 0 sin 8 — 6.183 cos? 0 sin 0 + 287.12 sin? 0 cos 6 


—233.45 sinf 0 cos 0 — 1531.8 sin? 6 cos 0 + 820.4 sinf 6 cos 6 + 3553.6 sin” 0 cos 0 
—734.76 sin? 0 cos 6 — 3831 sin? 8 cos 6 + 205.48 sin! @ cos 0 + 1546.8 sin!! 6 cos @ 
—220.48 sin!” 6 cos 0 + 137.34 sin" 6 cos @, 


ARCion(@L) — 7 389 cos 0 + 7.732 cos 0 sin 0 — 3.117 cos @ sin 8 — 128.12 sin? 8 cos 6 
—113.25 sin? 6 cos 0 + 515.16 sin? 0 cos @ + 374.15 sinf 0 cos 0 — 875.2 sin’ 0 cos 0 
—327.6 sin’ 6 cos 0 + 603.8 sin? 6 cos 0 + 83.061 sin!” 6 cos @ — 132.48 sin”! 6 cos 0 
—88.426 sin? 6 cos 0 + 65.01 sin'* 0 cos 0. 


Then, their changing rates are shown in Figures 7 and 8. 


ƏJRC( 4 LƏL 


AIRC 4 Lyla 0 
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Figure 8. The changing rates of the anisotropy. 


7. Discussion 


In this study, we investigated the anisotropy and scale effect of indeterminate JRC values by 


some unary NN functions and two-variable NN function, respectively, based on an actual case 


analysis. Generally speaking, the JRC characteristics reflected by the two-variable NN function 
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are more comprehension and offer more details, effectively indicating the anisotropy and scale effect 
of indeterminate JRC values simultaneously, although the NN function itself is more complex. 

Obviously, the proposed NN functions can not only express the JRC values with indeterminate 
information effectively and reasonably instead of the crisp values of JRC in the existing 
literature, but also indicates the anisotropy and scale effect of the indeterminate JRC values 
completely /simultaneously. However, the existing classical analysis methods cannot express and deal 
with such types of indeterminate JRC values and may lose much more information. Then, the main 
advantage of NN functions is that the NN functions of JRC can express and deal with the indeterminate 
JRC values instead of the crisp JRC values and contain much more useful information under 
indeterminate (imprecise, unsure, and even completely unknown) environments. Hence, we use 
a NN function instead of any classical function/crisp function of JRC to express the indeterminate 
function of JRC because any NN function can approximate the average value function/crisp function 
of JRC. When there is no indeterminacy in the NN function, the NN function is degenerated to 
the crisp function. Since, in the real world, we cannot always compute or provide exact/crisp values 
for the JRC characteristics, we can at least approximate them by NN functions. However, any NN 
function can contain and approximate the average value function/crisp function in existing classical 
methods. Obviously, the NN function provides one effective way for passing from the classical 
crisp functions to the NN functions in analyzing the anisotropy and scale effect of JRC under 
indeterminate environments. 


8. Conclusions 


Based on the concept and operations of NNs, the two NN functions JRC(@) and JRC(L) 
were developed for the first time to express the anisotropy and the scale effect of indeterminate 
JRC values, respectively. To indicate the anisotropic characteristic of JRC, the NN parameter 
was defined to quantify the anisotropic degree. Additionally, the derivative functions JRC’(@) and 
JRC’(L) were utilized to describe their changing rates in various measurement orientations and sample 
sizes. Finally, a two-variables NN function JRC(@, L) was presented to express the indeterminate 
information of JRC values in the sample sizes and measurement orientations and to investigate 
the changing rates by its partial derivative functions. It is obvious that the NN functions can effectively 
express the uncertainties of JRC and provide a new effective way for the JRC expression and analysis 
under indeterminate environments. In future work, we have to extend the NN function to the shear 
strength expression and conduct analyses under indeterminate environments. 
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